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Abstract. Using molecular dynamics simulations we investigate the thermodynamic of particles interacting 
with a continuous and a discrete versions of a core-softened (CS) intermolecular potential composed by 
a repulsive shoulder. Dynamic and structural properties are also analyzed by the simulations. We show 
that in the continuous version of the CS potential the density at constant pressure has a maximum for a 
certain temperature. Similarly the diffusion constant, D, at a constant temperature has a maximum at a 
density p max and a minimum at a density p m in < pmax, and structural properties are also anomalous. For 
the discrete CS potential none of these anomalies are observed. The absence of anomalies in the discrete 
case and its presence in the continuous CS potential are discussed in the framework of the excess entropy. 

PACS. XX.XX.XX No PACS code given 

1 Introduction However, experimental results have shown that for water 

at temperatures approximately below 10°C the diffusion 



Water is an anomalous substance in many respects. While 
most liquids contract upon cooling, for water the specific 
volume at ambient pressure starts to increase when cooled 
below T — 4°C at atmospheric pressure pQ. This effect 
is called density anomaly. Besides the density anomaly, 
there are more than sixty other anomalies known for wa- 
ter [3J. The diffusivity is one of them. For normal liquids 
the diffusion coefficient, D, decreases under compression. 



coefficient increases under compression and has a maxi- 
mum [3]. The temperature of maximum density (TMD) 
line inside which the density anomaly occurs, and the line 
of maximum in diffusivity are located in the same region 
of the pressure-temperature (P-T) phase diagram of wa- 
ter. Simulations for water also show thermodynamic and 
dynamic anomalies. The simple point charged/extended 
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(SPC/E) model for water exhibits in the P-T phase di- 
agram a TMD line. The diffusion coefficient has a max- 
imum and a minimum that define two lines at the P-T 
phase diagram, the lines of maximum and minimum in 
the diffusivity coefficient 4, 5 , 6J . Similarly to the experi- 
mental results, the TMD and the lines of maximum and 
minimum in the diffusion are located at the same region at 
the P-T phase diagram for the SPC/E model. Errington 
and Debenedetti [5] and Netz et al. [4] found, in SPC/E 
water, that there exists a hierarchy between the density 
and diffusion anomalies as follows. The diffusion anomaly 
region, inside which the mobility of particles grow as the 
density is increased, englobes the density anomaly region, 
inside which the system expands upon cooling at constant 
pressure. This observation is supported by experimental 
results [3]. 

Realistic simulations of water p3[Hl[S] have achieved 
a good accuracy in describing the thermodynamic and 
dynamic anomalies of water. However, due to the high 
number of microscopic details taken into account in these 
models, it becomes difficult to discriminate what is es- 
sential to explain the anomalies. On the other extreme, a 
number of isotropic models were proposed as the simplest 
framework to understand the physics of the liquid state 
anomalies. From the desire of constructing a simple two- 
body isotropic potential capable of describing the compli- 
cated behavior present in water-like molecules, a number 
of models in which single component systems of particles 
interact via core-softened (CS) potentials have been pro- 
posed. They possess a repulsive core that exhibits a region 
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of softening where the slope changes dramatically. This 
region can be a shoulder or a ramp [l0l[Tl"l[T2llf 3ill4ifT5l 
[THltl 7 , 18, 19, 2 OJ . Ramp and continuous shoulder- like po- 
tentials exhibit thermodynamic, dynamic, and structural 
anomalies 21,19,20 . However the square discontinuous 
shoulder shows no thermodynamic anomaly in three di- 
mensions |llj . All these potentials have in common the 
presence of two representative repulsive scales in the po- 
tential, a and ax, where the closest scale, a < a\, has the 
higher potential energy, U(a) > U{a{). 

One question that arises in this context is why a con- 
tinuous shoulder potential like the one described by de 
Oliveira et al. [TW2"U] has the anomalies while the square 
shoulder described by Franzese et. al. [11] has no anomaly? 
In order to shade some light in the reasons for the pres- 
ence of density anomalies in CS potentials, both the dis- 
continuous shoulder potential and the continuous version 
are analyzed in the framework of the excess-entropy-based 
formalism recently introduced by Errington et al. |22J. 
Within this approach the presence of the density anomaly 
is related to the density dependence of the excess entropy, 
Sox- We will follow this surmise and compute the excess 
entropy for both the discontinuous and the continuous 
shoulder potentials. From the analysis of s cx for the dis- 
continuous model we are able to derive a simple argument 
for the presence or absence of anomalies in CS potentials. 

The outline of the paper is as follows. We present the 
details of the two models in Sec. [H details of the simu- 
lations, and the P-T phase-diagram for both models are 
presented in Sec. [H In Sec. H] comparisons between the 
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behavior of the excess entropy of the continuous and dis- 
continuous models are made and from that a condition for 
the presence or absence of anomalies in CS potentials is 
proposed. Conclusions end this sections. 

2 The models 

The first model we study here consists of a system of N 
particles of diameter cr, inside a cubic box whose volume is 
V, resulting in a density number p = N/V [19j interacting 
with a continuous shoulder potential given by 
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where U*(r) — U(r)/e. The first term of Eq. {T) is a 
Lcnnard- Jones potential of well depth e and the second 
term is a Gaussian centered on radius r = r$ with height 
a and width c. In a previous publication we have studied 
this model with setting a = 5, tq/u = 0.7 and c = 1 (see 
Fig- H) [IH]- This potential has two natural length scales: 
one close to the hard core, cr, and another at a further 
distance where the potential has its lower value. This last 
length we call <7\ . 

The second model we study here is a system of N 
particles of diameter cr, inside a cubic box whose volume is 
V, resulting in a density number p — N/V but interacting 

with a discontinuous shoulder potential given by 

oo r < a 

U*(r) = { 1 di > r > cr , (2) 
r > a\ 



U (r) 




Fig. 1. Interaction potential from eq. (JTJ with parameters a - 
5, ro/a — 0.7 and c = 1, in reduced units. 
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Fig. 2. Interaction potential from eq. ((2} with parameter 
d/cr = 1.75 in reduced units. 

where U*(r) = U(r)/e. This potential has two natural 
length scales: the hard core distance, cr, and the outer 
core, o\. Here we analyze the case u\ja = 1.75 illustrated 
in Fig. [21 Here we use dimensionless pressure, P* , tem- 
perature, T*, and density, p*, that are given in units of 
cr 3 /e, fce/e, and cr 3 , respectively and ks stands for the 
Boltzmann constant. 

3 Details of simulations 

For the continuous shoulder potential we performed molec- 
ular dynamics simulations. Details of the simulation can 
be found in ref [TH]- Figure ([3]) shows the P-T phase- 
diagram we have obtained in a previous publication |19j . 
The isochores have minima that define the temperature 
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Fig. 3. Pressure-temperature phase-diagram obtained for the 
continuous shoulder potential. From bottom to top p* = 0.100, 
0.103, 0.105, 0.107 0.110, 0.113, 0.115, 0.117, 0.120, 0.123, 
0.125, 0.127, 0.130, 0.132, 0.134, 0.136, 0.138, 0.140, 0.142, and 
0.144 are shown. The solid line illustrates the TMD and the 
dashed lines show the boundary of the diffusivity extrema. 




Fig. 4. Pressure-temperature phase-diagram for the discontin- 
uous shoulder potential. From bottom to top p* = 0.08, 0.09, 
0.10, 0.11, 0.12, 0.13, 0.14, 0.15, 0.16, 0.17, 0.18, 0.19, 0.20, 
0.21, 0.22, 0.23, 0.24, 0.25, and 0.26 are shown. 

of maximum density. The TMD line encloses the region 
of density (and entropy) anomaly. We also have studied 
the mobility associated with the potential described in 
Eq. ([T]) [TH]. The diffusion was also calculated using the 
the mean-square displacement averaged over different ini- 
tial times. The behavior of D as a function of p* goes as 
follows. At low temperatures, the behavior is similar to 
the behavior found in SPC/E supercooled water 0j. The 
diffusivity increases as the density is lowered, reaches a 
maximum at pDm&x (and PD max ) and decreases until it 



reaches a minimum at pDmin (and Pomm)- The region in 
the P-T plane where there is an anomalous behavior in 
the diffusion is bounded by (T Dmin , Pomm) and {T Dms , x , 
PDmax) and their location is shown in Fig. (J3J). The re- 
gion of diffusion anomalies (T Dmax , Pomax) and (T Dmin , 
Pomin) lies outside the region of density anomalies like in 
SPC/E water g]. 

In order to simulate the discrete potential showed in 
Fig. ([2]) we used the collision driven molecular dynamics 
techniques [23) . 500 particles were put into a simulation 
box with periodic boundary conditions and the rescaling 
velocities scheme was used for every 1000 collisions un- 
til equilibration time to achieve the desired temperature. 
After thermalization particles were allowed to move under 
microcanonical ensemble. In the collision driven molecular 
dynamics simulation, kinetic energy has to be rigorously 
conserved. Hence, no special mechanism is necessary in 
order to simulate a desired temperature and the NVE en- 
semble becomes the natural choice. The equilibration and 
production times in reduced units were 350 and 650 re- 
spectively. Figure ([!]) shows the P-T phase-diagram for 
the discontinuous shoulder potential. The isochores for 
different temperature and pressures show no minima so 
no density anomaly is present. 

4 Excess entropy and anomalies 

Why the discontinuous shoulder potential has no water- 
like anomalies and its continuous counterpart has all of 
them? We can gain some understanding about that by 
analyzing the density dependence of the excess entropy. 
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Errington et al. have shown that the density anomaly is 
given by the condition £ cx = (ds ex /<91n p) T > 1 [22]. Here 
Sox is the excess entropy and is approximated by its two- 
body contribution, 

,S2 = —2np J [g(r) In g(r) — g(r) + 1] r 2 dr . 

The radial distribution function, g(r), is proportional 
to the probability to find a particle at a distance r to 
another particle placed at the origin. Errington et al. [22j 
have also suggested that the diffusion anomaly can be pre- 
dicted by using the empirical Rosenfeld's parameterization 
They found the condition £ 2 > 0.42 for a diffusion 
anomalous behaviour. They also claim that £% > is a 
good estimative for determining the region where struc- 
tural anomaly occurs |22j . 

In order to understand the differences between the con- 
tinuous and the discontinuous shoulder potentials we test 
the excess entropy criteria described above in both poten- 
tials. The radial distribution functions for different tem- 
peratures and densities for both potentials were then ob- 
tained by the molecular dynamic simulation method de- 
scribed in Sec. [3 

Figure (O illustrates the two-body contribution of ex- 
cess entropy for the continuous potential given by Eq. ([T]) . 
S2 is negative and its slope changes from positive to neg- 
ative what indicates the presence of structural anomaly. 
Figure ([6]) shows the behavior of £2 with density for a 
fixed temperature for the continuous model. The horizon- 
tal lines at £\ — 0, 0.42, and 1 indicate the threshold 
beyond which there are structural, diffusion, and density 
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anomalies respectively. In accordance with Fig. §5§ the 
density anomalous region shown in Fig. ([6]) occurs in an 
interval of density smaller than the interval where the dif- 
fusion is anomalous. The two-body excess entropy also 
show the presence of structural anomaly what corrobo- 
rates results of simulations of structural parameters [5D] . 

Figure ([7]) illustrates the two-body excess entropy of 
the discontinuous potential given by Eq. @ . S2 is negative 
but its slope is for almost all temperatures and densities 
negative. The filled circles show the region where the sys- 
tem crystallization occur. Figure (JSj) shows the behavior of 
£2 with density for fixed temperatures for the discontinu- 
ous model. The line £\ = 1 is never crossed so no density 
anomaly should be expected what is in good agreement 
with the Fig. 0$. The line £\ = 0.42 is also never crossed 
what indicates that diffusion anomaly is not expected in 
the discontinuous model. This is also in agreement with re- 
sults for similar step potentials where no diffusion anomaly 
is found for large steps [25] . The line £2 — is crossed 
for temperatures T* = 0.17 and 0.2 what would suggest 
the presence of structural anomaly. However this has to 
be taken with a grain of salt since Errington et al. [22] 
demonstrated that S2 overestimates the region of anoma- 
lies. In this sense, a detailed study of translational and 
orientational order parameters [20 is necessary prior to 
any affirmative on the presence of structural anomaly for 
this discontinuous shoulder. 

Even thought the excess entropy criteria is quite use- 
ful for predicting if the anomalies would be present for a 
certain potential it does not provide an easy and intuitive 
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method for understanding why the continuous shoulder 
has anomalies and the discontinuous one does not have. 
In principle, both potential exhibit similar characteris- 
tics. Both potentials have two repulsive scales and conse- 
quently the radial distribution function in both cases has 
two peaks, one close to the hard core and another close to Fig- 5- Pair contr ibution of excess entropy, s 2 , against density 
the distance r = u x as illustrated in Fig. © for the con- f or T* = 0.25, 0.30, 0.35, 0.40, 0.45, 0.50, 0.70, 1.00, and 3.00 
tinuous potential and in Fig. (JTUJ) for the discontinuous from bottom to top for the continuous shoulder model (Fig. iJTJ 
potential. A closer look at the radial distribution reveals . Circles are simulated data and lines are fifth order polynomial 
an important difference between the two cases. For the fit from data. 



continuous potential for densities and temperatures in the 
region where anomalies occur g(r) grows with density for 
r/c « 1 and decreases with increasing density for r rs <7i 
[see Fig. (j9]) for example] . For the discontinuous potential 
the g(r) increases with density both at the hard core and 
at r = o\ for any temperature and density [see Fig. (|10[) 
for example]. Notice that the radial distribution function 
in both cases has significative changes with the change in 
density at the two natural scales. 



uously. As it was shown by Netz et al. [25 for the steps 
potential, anomalies would only be observed if the discon- 
tinuity in U* would be below a certain threshold. 

Now we shall test if this simple hypothesis is in agree- 
ment with Errington et al. criteria. First we compute S2 



So = 



ds 



dlllpy rp 

s 2 - 2-Kp 2 I In g(r) 



2 / ' -^ (r W 



dp 



Finally we propose that a two scale potential has anoma- In this expression, the first term is always negative for all 

lies for some temperature and densities if dg{r)/dp > two scale potentials. The integral in the second term, as 

for r « a and dg(r)/dp < for r ss o\. If this require- we have observed above, is dominated by the the values 

ments would not be fulfill for any temperature and density at the two scales, a and o\. 

no anomaly would be present. This proposition is based For the continuous potential in the region where the 

in the physical picture that within the anomalous region anomaly is present In g(r « a) < 1 while In g(r ~ a±) > 1. 

for a fixed temperature an increase is density implies an Also dg(a)/dp > while dg(o\)/dp < 0. Consequently 

increase in the number of particles close to the hard core, the second term in Eq. ([3]) is positive. This allows for a 

This particles move from the distance a\ to a. For a dis- zero or positive values of £2 for appropriated densities 

continuous potential this requires an activation energy of and temperatures. Therefore our criteria is in accord with 

U = e while for the continuous case it can be done contin- Errington et al. criteria. 
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Fig. 6. £2 = (dsz/dXn p) T versus density. Following the iso- 
chore p* = 0.10, from top to bottom, temperatures T* = 0.25, 
0.30, 0.35, 0.40, 0.45, 0.50, 0.70, 1.00, and 3.00 are shown for 
the continuous shoulder potential. 




Fig. 7. Pair contribution of excess entropy, S2, against density 
for T* = 0.15, 0.17, 0.20, 0.23, 0.25, 0.30, 0.35, 0.40, 0.45, and 
0.50 from bottom to top for the discontinuous shoulder model 
(Fig.©. Open circles are simulated data and lines are fifth 
order polynomial fit from data. Solid circles show the limit of 
crystalization. 
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Fig. 8. £2 = (pS2/ 'din p) T versus density for the discontinu- 
ous shoulder potential. Following the isochore p* — 0.13, from 
bottom to top, T* = 0.15, 0.17, 0.20, 0.23, 0.25, 0.30, 0.35, 
0.40, 0.45, and 0.5 are shown. 
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Fig. 9. Radial distribution function for the potential given by 
Eq. JTJ for T* = 0.25 and densities p* = 0.12, 0.14, 0.16, and 
0.18. Arrows indicate the direction of increasing p* . 
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Fig. 10. Radial distribution function for the potential given 
by Eq. (2j for T* = 0.25 and densities p* = 0.12, 0.14, 0.16, 
and 0.18. Arrows indicate the direction of increasing p* . 

In resume, in this paper we have calculated the ex- 
cess entropy and its derivative for both continuous and 
discontinuous shoulder potential. For the continuous case, 
using the Errington et al. criteria indicates that this po- 
tential has density, diffusion and structural anomalies as 
we have shown in previous publications. For the discon- 
tinuous potential the criteria indicates that no thermo- 
dynamic and dynamic anomalies are present and its not 
conclusive for structural anomalies. Direct calculations of 
the P-T temperature phase-diagram confirms the excess 
of entropy prediction. On basis of these results we pro- 
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pose a criteria for predicting if a two scale potential has 
or not anomalies. Our criteria provides a simple picture 
for the anomalies not being observed in the discontinuous 
shoulder potential. 
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